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Nomenclature 


mass matrix for substructure In {x) physical coordinates 
stiffness matrix for substructure in {x} physical coordinates 
matrix of substructure eigenvectors In {x> coordinates 
mass matrix for total structure In {x} physical coordinates 
stiffness matrix for total structure In {x} physical coordinates 
matrix of eigenvectors for total structure in {x} physical coordinates 
generalized mass matrix for total structure In {q} modal coordinates 
generalized stiffness matrix for total structure In {o} modal coordinates 
matrix of eigenvectors for total structure In {q} modal coordinates 
dynamic transfonnatlon matrix 

transformation matrix defining the relationship of the reduced 
coordinates, {qR> , to the kept coordinates {qK}. 

physical coordinates 

generalized modal coordinates 

system circular frequency (rad/sec) 

substructure circular frequency (rad/sec) 

reduction circular frequency (rad/sec) 

number of degrees of freedom in structure 

Subscripts 

refers to the 1^*^ subset, term, or substructure 
Incremental mass from coupling spring 
Incremental stiffness from coupling spring 
kept coordinates 
reduced coordinates 


11 



Superscripts 


kept coordinates 
reduced coordinates 
attachment coordinates 

interior coordinates not attached to any other substructure 

revised value 

second time derivative 

particular submatrix partition 


I 


iii 



SECTION 1 
INTRODUCTION 

This document describes a Fortran IV computer program used to perform modal 
synthesis of structures by stiffness coupling using the dynamic transformation 
inethod. The program has been named SCAMP (Stiffness Coupling ^proach Modal- 
Synthesis Program). The program begins with the entry of a substructure's 
physical mode shapes and eigenvalues or a substructure's mass and stiffness 
matrix. If the mass and stiffness matrices (100 degrees of freedom [DOF] 

Limit) are entered the eigen problem for the individual substructure Is 
solved. Provisions are Included for a maximum of 20 substructures which may 
be coupled by stiffness matrix springs (90 DOF/spring). Each substructure 
can have a number of OOF, except that for OOF greater than 100, vector sets 
having maximum row and column size of 100 have to be generated prior to 
entering SCAMP. The substructures .are then coupled together via coupling 
springs, and the dynamic transformation is used to reduce the size of the 
eigen problem. The total number of modes treated by the program is 300 
consisting of 100 kept coordinates (maximum eigen value size) and 200 
coordinates reduced by the dynamic transformation. For user flexibility, 
five major entry points have been included with restarting capabilities in 
SCAMP. Tiering can also be used with SCAMP. The assembled solution can be 
identified as a substructure and added to other substructures if desired. 

This procedure can be repeated indefinitely but care has to be exercised on 
the conditioning of the system mode shapes and frequencies. 

Input data for SCAMP Is accomplished mainly by the "READ" FORMA subroutine. 
Output data to be saved for re-starting SCAMP is written by special subroutines. 
Files containing data should be copied to an output data tape ($}. Modes and 
frequencies for the coupled structure are output on a user specified file (or 
tape) and can be read using the "READ" subroutine. 
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SECTION 2 


THEORETICAL DISCUSSION 
2.1 BASIC THEORY FOR STIFFNESS COUPLING 

The stiffness coupling method of modal synthesis assembles the complete structure 
in the same manner as the displacement method for structural analysis. The total 
structure may be represented by a number of substructures connected through flexible 
links. Each substructure is analyzed without the flexible links to determine the 
component vibration modes with free attachment coordinates. The flexible links are 
represented by a stiffness matrix relating the interface forces frwi one set of 
substructure attachment coordinates to another. 

The method of substructuring for stiffness coupling may best be illustrated by 
considering a total structure consisting of only two substructures. The general 
undamped equation of motion for the ith substructure in terms of its 
mass matrix, [ ^i], and stiffness matrix, [JS’Ji is given by 

-t- -O (1) 

in which the coordinates y describe physical motions of the mass points. Each 
substructure has two sets of coordinates which will be referred to as attachment 
coordinates, { }, and internal coordinates, { >7 The attachment coordinates are 

those degrees of freedom (DOF) which are connected to another substructure via a 
stiffness matrix or coupling spring. The internal coordinates are those DOF which are 
not connected to any other substructure. It is important to note that for stiffness 
coupling the coordinates belonging to one particular substructure are not common to 
any other. If n* represents the size of the ith substructure and m the size of 
the total structure comprised of f substructures, then n will be given by 
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Having defined two sets of coordinates for each substructure, Eo. (1) may be written 
in partitioned form as 
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Now consider the total structure to be described by a mass matrix, [ftfl, and a 
stiffness matr*A [K-j-], such that 




where 


[Mr] 


II I 

My ( 0 






I 


( 5 ) 


2-2 



If we describe the connect<ng structure between (x^) and (x*) by a mass matrix, 
[m^3, where 
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and a stiffness matrix relating the nodal force, {F*}, to (x*) by the 

equation 
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Then [M^] and [K^] may be written In partitioned form using the submatrices defined 
by Eqs, (3), (6) and (7): 
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and recognizing the partitions belong to each of the substructures, Eo. (4) may 
be written as 
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For each substructure defined by Eq. (1), a set of eigenvalues, 

and a set of mass normalized eigenvectors, [^^], can be obtained such that 
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Using the results from Eq. (11), we can now express Eq. (10) in terms of a 
set of generalized modal coordinates, {q}. The coordinate transformation is 
given by 


where 


h«ri w* I 

I'i'fivj 

*'•1 n<i n*, 

W . Kflc] 


V\»H 

Substituting the transformation for (x) into Eq. (10) and premultiplying by 
yields 
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Solution of Eq. (14) will result in a set of eigenvalues* and their cor- 

responding eigenvectors, [If]. Because of the similarity transformation used to 
obtain Eq. (14) from Eq. (10), the eigenvalues of Eq. (10; will be equal to those 
of Eq. (14) and the corresponding eigenvectors, Eq* (^0) will be given by 

[ ^*1 * [*] ( 15 ) 

Eq. (14) represents’ the most general form of the equation of motion for stiffness 
coupling. This equation is g^erally solved by partitioning the {q} coordinates 
into two groups, kept and truncated. The truncated coordinates correspond to the 
high frequency substructure modes and are completely omitted from the equation of 
motion. Those degrees of freedom remaining, the partitioned set of kept coordinates, 
determine the final reduced size of the eigenvalue problem to be solved. 

The general form of Eq. (14) may be simplified further by including the correct 

m . partitions from Eq. (6) in each corresponding [m.j] at the substructure 
^ 1 

level. This is reasonable because we are assuming that there is no inertial coupling 
between substructures. This will result in [M.] » 0 ard Eq. (14) reduces to 
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After solving for the w^'s and ^^'s, the only lengthy calculation left to be per- 
formed In order to obtain Eq. (16) Is the matrix tripie-product involving [Kqpl]- 
If we partition each [(^^] row-wise in terms of its n^ and n| coordinates such that 


Hi ‘Hi 

Eq. (13b) can be written as 
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By forming the triple product 'Re- 

factoring the result in terms of the matrix partitions, the resultant form of the 
triple product may be expressed as 
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and the final form of [K] in Eq. (17) may be expressed as 
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2.2 DYNANiC TRANSFORMATION 

As a result of omitting the higher substructure modes, the solutions from the 
truncated Eq. (16) will have errors introduced. The truncation errors can be 
greatly diminished by including the modes that would have been truncated through a 
dynamic transformation. Instead of truncating or omitting modes, all modes can be 
included through a transformation that relates the "reduced" modes not contained 
explicitly In the solution to the mcxies that are "kept." If n? corresponds to an 
exact eigenvalue of Eq. (16), the relationship between the eigenvalue and its eigen- 
vector may be expressed in terms of the kept, {q*^}, and reduced, {q*^}, coordinates as: 



K«l 

where the {q*^} corresponds to those modes previously truncated. If we designate nj^ 
as the total number of modes kept from all the substructures and n^ as the total 
number of modes reduced, then 

n = n|( + n^ (23) 

7 2 

Expanding Eq, (22) into two equations for some general frequency, p = a-j , yeilds 
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Solving Eq. (24b) for {q*^} in terms of {q*^} gives 

h'i- W!»‘i 

hy«l 

where 

Qr] , - Ck"- p‘I]“ [<''*] 

Using Eq. (25) for some '‘reduction'-* frequency, we can write 

b-]m 

h&l h«» *^1it** 

The dynamic transformation matrix, [T], is then defined as 
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The reduced equation of motion Is obtained directly by substituting the coordinate 
transformation 

f\K\ 


into Eq. (16) and pre-multi plying by the transpose of [T]. The reduced generalized 
mass and stiffness matrices can be written into the partitioned forms given by 
Eqs. (22) and (28) • 


CM*] = [I] [R]^[R] 


!1k"]= [K**Jf2[/<*-JW+[Rr[K-j[R] 
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Conventional methods of determining eigenvalues may be applied to the reduced 
equation of motion to obtain a set of eigenvalues, and a corresponding 

set of mass normalized eigenvectors, From the coordinate relationship 

defined by Eq. (25), the reduced eigenvectors, corresponding to the {q>"} 

reduced coordinates are given by 


and the physical eigenvectors for the total solution will 


where 





be given by 

(32) 

(33) 


This solution will be exact for any ^ which is the same as the reduction frequency, 
p , used in developing [T]. 



SECTION 3 


PROGRAM ENTRY POINTS 

In order to provide user flexibility and at the saire time ir.inimize computation of basic 
data changes, five major entry points have been established for SCAMP: 

1. Basic substructure data entered, subsysten eigensolutions. 

2 . Coupling spring stiffness data entered and stiffness contributions 
calculated. 

3. Selection of modes to be kept/reduced/ truncated, generalized mass and 
stiffness matrices calculated. 

4. p2 value for dynamic transformation entered, eigensolution for system 
using the dynamic transformation. 

5. Calculation of physical eigenvectors for total systen, and user "print" 
options. 

Entry into the program is accomplished by designing a specific entry point (EP). 
Termination of theprogram is accomplished by designating the last EP the user 
desires to execute. 

The Figure 3-1 Flow Chart is included to show the general flow of the program. 

Entry Point 1 

The mass and stiffness matrices [iHj], [k^3, or mode shapes and eigenvalues [1^]. 

CX-j] are read into SCAMP at EP-1 . The substructures are defined in the program by 
a user supplied number which ranges from 1 to 20. Since each substructure is to be 
Identified in the program by a distinct number, their input may be in any order. 

Each fm^], [k. ] is limited to IOC DOF. If modes and frequencies are input for 
substructure definitions, there is no limit on the physical degrees of freedom for 
the substructure. However, if more than 100 physical DOF are used, the matrix of 
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Program Flow Chart 







Entry Point 3 

At this EP, the user specifies which modes are to be kept and reduced. The input 
matrix IV defines how many modes are kept and reduced for all of the substructures. 

The elements are read in groups of three where the first element corresponds to 
some substructure identification number "ISUB", the second element the number of 
modes to be “Kept" for ISUB, and third element the number of modes to be "reduced" 
for ISUB. All other modes will be truncated. The restrictions for the total kept 
and reduced modes are given by 

maxsub 

I KEEP (1. i) < 100 
1 * 1 
maxsub 

I KEEP (2, i) <200 
i * 1 

The generalized stiffness matrix is then assembled in its kept and reduced partitions. 
Only the ISUB's specified by IV will be coupled together even though there may be 
several more substructures and springs residing in the basic data set. 

Entry Point 4 

The reduction frequency, p, is entered at EP-4 (where p2 = LAMDA0), and the dynamic 
transformation is aoDlied to obtain the reduced mass and stiffness matrices. The 
eiaenvalues for the couoled substructures are now obtained. 

Entry Point 5 

This entry point calculates the coupled physical eigenvectors which are printed in 
substructure groups. They are printed in ascending order from the originally assigned 
substructure number. The substructure partitions, upon output, are resequenced from 
one to some number "m". This facilitates re-starting SCAMP with the coupled system as a 
substructure having "m" partitioned physical degrees of freedom. The user must 
specify how many modes are to be printed by specifying the first and last mode number 
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of the group of modes desired. Considerable computer time may be saved by printing 
a few of the modes. The user may desire to re-enter the program to select more 
modes. Also participation factors can be printed out by specifying the first and 
last mode number of the modes desired. Minimum values for these participation 
factors can be specified. Absolute values less than this will be set equal to zero 
and not printed. An output tape option exists which specifies the number of modes 
and frequencies to be put on an output user tape. The tape will contain coupled 
eigenvalues, frequencies in Rad/Sec, frequencies in Hz, and eigenvectors and can be 
read by the "read" subroutine. 


a 
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SECTION 4 


INPUT/OUTPUT DATA 


4,1 INPUT DATA 

This section describes the necessary Input data required to execute SCAMP and Is 
ordered by entry points. The data cards and x.he variables appearlny on them are 
listed In sequential order with the Input format or FORMA subroutine specified. The 
definitions for the variables specified on the cards are given to clarify their 
meaning. Variables followed by A.V. (Assigned Value) specifies that this particular 
variable has been pre*ass1gned a specific given value. If the A.V. specifies the 
desired user option, the user may delete defining the variable on the input list. 

Explanation for “READ" formatted data given at end of Input data section. 


Input Data for Starting Program - Entry Point 0 
Card Number Input Data 

1 IRUN, UNAME, Nl. T1 

2 TITLEl 

3 TITLE2 

4 lENTR, lEXIT, IFIRST, F|D 


Fornat 

(A6. 4X, 3A6, 2X, 12, A6) 
(12A6) 

(12A6) 

$ ENTR 


Definitions 


IRUN 


* Run Number 


UNAME 

Nl 


T1 

TITLEl 


User's Name 

Write-Tape-No. (Eg, 21). If coupled modes and eigenvalues are to 
be output on an uninitialized tape. The Write-Tape No. must be 
specified here to Initialize the tape. Leave blank for no tape 
Initial Izatlon. 

Write-Tape-ID(E6, T12345). Leave blank if Nl not specified. 

First Title 


TITLE2 


Second Title 


$ ENTR 


* Namelist for Card 4 Input 
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.. ... I™. i .......... .. l.. . . I - I I .....4 ... i 


lENTR * Entry Point (EP) for entering SCAMP. There are 5 EP's in all defined 

by the integers from 1 to 5. 

lEXIT * Exit Point for terminating SCAMP. The exit points are identical to 

the entry points with program termination occurring after the 
execution of the EP specified by lEXIT. The program returns to 
Card 1 above. If IRUN equals the literal "STOP", the program 
termi nates. 


IFIRST 


0 (A.V.)*, denotes that SCAMP has been executed before and a restart 
tape containing previously generated data will be used. 


m 

1, required for entering a new problem 

into SCAMP 

for the first time. 

F(|D 

l.E-16 (A.V.), final off-diagonal value 
matrix using the method of Jacobi. 

for diagonalizing a syiranetrlc 

*The Assigned Value (A.V.) given by SCAMP. If the A.V. 
the usar may delete the variable from the input list. 

specifies 

the correct option. 

Entry Point 

2 



Card Number 

Input Data 

Format 


1 

ISUB. NMteS, KEPM0D(2), MASSD, 
MASSWT, lALTER, M$DES, IRIGIO, 
IR0T, NSUBS. lEND 

SENTRl 


*2 

DIRC$S (3 X 3) 

READ 


*3a 

MASS (NXN) or (IXN) 

READ 


♦4a 

STIFF (NXN) 

READ 


♦3b 

LAMBDA (IXM) 

READ 


*4b 

PHY [NSUBS*{NXM)] 

READ 


N.B. Card 1 

with lEND * 1 is the last card for substructure input. 



Definitions 



SENTRl = 

NAMELIST for Card 1 Input. 




ISUB = Substructure identification number which ranges from 1 to 20, Each 

substructure must be identified by a different number. 

NM^DES * The number of degrees of freedom for the mass and stiffness matrices 

of substructure "ISUB*'. If physical modes are to be read, NMODES 
equals the number of modes of "ISUB". (NMODES^j^j^ = 100) 
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i 1 1 I I = I ■ f 


KEP^^D(1) 

KEPr*^D(2) 

MASSD 

MASSWT 

lALTER 

M^DES 

IRIGID 

IR0T 

NSUBS 

lEND 

O1RC0S 

MASS 


■ NM0DES (A.V.), number of modes to be used by SCAMP for substructure 
representation. Values less than NMODES v/111 truncate the hichest 
NM0DES-KEPMOD(1) substructure's modes. 

■ NM0DES (A.V.), nuiri}er of modes to be printed. 

* -1, deletes printing of substructure modes. 

■ 1 (A.V.), diagonal mass matrix option where the diagonal mass is a 
1XNM0OES vector. 

* 0, denotes generalized mass matrix option (NF^DES x NMfOES) 

* 1 (A.V.), weight matrix to be read Instead of mass 

■ 0, mass matrix option 

* 0, (A.V.), alter option. The current "ISUB" Is a new substructure 
being added to the basic data set. 

» 1, the current "ISUB" has been read into SCAf^P previously. This option 

will replace the "old" ISUB with new data. 

* 0 (A.V.), option for entering substructure mass and stiffness matrices 

* 1, option for entering substructure eigenvalues and eigenvectors. 
Options referring to mass and stiffness data are Ignored. 

= 0 (A.V.), use eigenvalues as caluated or read Into SCAMP. 

» NR, zero-out first NR eigenvalues before saving on tape. 

* 0 (A.V.), do not rotate substructure 

* 1 , rotate substructure 

■ -1, rotate substructure and print un-rotated eigenvectors. KEPM0D(2) 

» -1 suppresses all eigenvector printing. 

s 1 (A.V.), number of eigenvector sets within a substructure to be 
read In under MODES®! option for substructure "ISUB". Maximum row 
and column size for any one vector set Is 100. There Is no limit for 
NSUBS. 

■ 0 (A.V.), optional end-of-data card. Program continually returns 

to Card 1 until user soecifies lEND * 1$. Cards 2 through 4 ore then 
skipped and EP-1 terminated. 

* Input direction cosines matrix (3X3). Required only when IR0T.NE.C. 
Total number of substructure DOF's must be divisible by 3 and ordered 
for rotation in order to exercise this option. 

» Input mass matrix for "ISUB" (NXN) or (IXN) 
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STIFF 

LAMBDA 

PHY 

Entry Point 2 
Card Numbtr 
1 
2 

3 

4 

SENTR2 

IJSUBSm » 

IJSUBS{2) * 

NACORDO) « 
NAC0RD{2) - 
lALTER 

m 

I END 

SENTR2A.B - 
NASUB 


Input stiffness matrix for "ISUB” (NXN) 
Input eigenvalues for "ISUB" (IXM) 

Input eigenvectors for "ISUB" (NXM) 
(Mmax • %X " 


Input Data 

Format 

IJSUBS (2), NAC0RD(2), lALTER. lEND 

SENTR2 

NASUB, NA (NC) 

SENTR2A 

NASUB, NA (NC) 

$ENTR2S 

KCPL (MXM) 

READ 


Definitions 

NAMELIST for Card 1 Input. 

First "ISUB" coupled by KCPL. DOF for this substructure will be 
located In the KCPL-i^ partition of KCPL. "ISUB" is the user speci- 
fied substructure number. 

Second "ISUB" coupled by KCPL. DOF for this substructure will be 
located in the KCPL 22 partition of KCPL. 

Total number of attachment DOF's in IJSUBS(l). 

Total number of attachment DOF's in IJSUBS(2). 

0 (A.V.), alter option. The current KCPL defined by IJSUBS Is a new 
coupling spring being added to the basic data set. 

1, the current KCPL has been read into SCAMP previously. This 
option will replace the "old" KCPL with new data. 

0 (A.V.), optional end-of-data card. Program continually returns to 
Card 1 until user specifies lEND ■ IS. Cards 2 through 4 are then 
skipped and EP-2 tenninated. 

NAMELIST Input for Cards 2 and 3. 

Number of coordinates (elements) specified in "NA" vector. 
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KEPMfDd) 

K£PM$0(2) 

MASSO 

MASSWT 

lALTER 

M$DES 

IRI6ID 

IRCIT 

NSUBS 

lEND 

DIRCfS 

MASS 


* NM0OES {A.V.)i nuft^er of modes to be used by SCAMP for substructure 
representation. Values less than NMODES will truncate the hichest 
NM^OES-KEPMOO(l) substructure's modes. 

■ NM|0ES (A.V.), number of modes to be printed. 

* -1, deletes printing of substructure modes. 

* 1 (A.V.), diagonal mass matrix option v^ere the diagonal mass Is a 
UNHIDES vector. 

■ 0, denotes generalized mass matrix option (N^^DES x NMfDES) 

■ 1 (A.V.), weight matrix to be read Instead of mass 

■ 0, mass matrix option 

* 0, (A.V.), alter option. The current "ISUB" 1s a new substructure 
being added to the basic data set. 

* 1 , the current "ISUB" has been read into SCAf^P previously. This option 
win replace the "old" ISUB with new data, 

* 0 (A.V.), option for entering substructure mass and stiffness matrices 

* 1, option for entering substructure eigenvalues and eigenvectors. 
Options referring to mass and stiffness data are ignored. 

» 0 (A.V.), use eigenvalues as caluated or read Into SCAMP. 

» NR, zero-out first NR eigenvalues before saving on tape. 

» 0 (A.V.), do not rotate substructure 

■ 1 , rotate substructure 

» -1, rotate substructure and print un-rotated eigenvectors. KEPMCDfS) 

= -1 suppresses all eigenvector printing, 

» 1 (A.V.), number of eigenvector sets within a substructure to be 

read in under M00£S*1 option for substructure "ISUB". Maximum row 
and column size for any one vector set is 100. There is no limit for 
NSUBS. 

« 0 (A.V.), optional end-of-data card, '■’rogram continually returns 

to Card 1 until user specifies lEND » 1$. Cards Z through 4 are then 
skipped and EP-1 terminated. 

* Input direction cosines matrix (3X3). Required only when IRlT.NE.C. 
Total nunter of substructure DOF's must be divisible by 3 and ordered 
for rotation in order to exercise this option. 

» Input mass m.atrix for "ISUB" (NXN) or (IXN) 
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I 


STIFF ■ Input stiffness rnatrix for "ISUB" (NXN) 

LAMBDA * Input eigenvalues for "ISUB" (IXM) 

PHY * Input eigenvectors for "ISUB" (NXM) 


Entry Point 2 


Card Number 
1 
2 

3 

4 


Input Data 

Format 

IJSUBS (2). NACfRD(2), lALTER, lEND 

SENTR2 

NASUB. NA (NC) 

SENTR2A 

NASUB, NA (NC) 

SENTR2B 

KCPL (MXM) 

READ 


Definitions 


SENTR2 * NAMELIST for Card 1 Input. 

IJSUBS(l) ^ First "ISUB " coupled by KCPL. DOF for this substructure will be 

located in the KCPL^I partition of KCPL. "ISUB" is the user speci- 
fied substructure number. 


IJSUBS(2) » Second "ISUB" coupled by KCPL. DOF for this substructure will be 
located in the KCPL 22 partition of KCPL. 

NACUPD(l) * Total number of attachment DOF's in IJSUBS(l). 

NAC0RD(2) * Total number of attachment DOF's in IJSUBS(2). 


lALTER 


0 (A.V.), alter option. The current KCPL defined by IJSUBS is a nev. 
coupling spring being added to the basic data set. 


» 1 , the current KCPL has been read into SCAMP previously. This 

option will replace the "old" KCPL with new data. 

lEND * 0 (A.V.), optional end-of-data card. Program continually returns to 

Card 1 until user specifies lEND « IS. Cards 2 through 4 are then 
skipped and EP-2 terminated. 


SEI 1 TR 2 A.B . NAMELIST Input for Cards 2 and 3. 

NASUB * Number of coordinates (elements) specified in "NA" vector. 
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NA « Input vector identifying the attachment coordinate IK)F's. $ENTR2A 

selects coordinates from substructure IJSUBS(l) and $£NTR2B selects 
coordinates from 1JSUBS(2). If NSUBS = 1, the first element NA(1) 
will be equal to 1 and the remaining elements will define all the 
attachment coordinates In the substructure. If NSUBS>1 , all the 
attachn«nt coordinates may not occur In one eigenvalue partition of 
the substructure. In this case NA(1) designates the eigenvalue 
partition number and the remaining elements define the attachment 
coordinates in this partition. When NA(1) is given a negative sign^ 
the namelist is repeated in order to define attachment coordinates 
from another partition. A positive sign is assigned to NA(1) for 
the partition that completes the specification of attachment 
coordinates. 

The part of the vector identifying the attachment coordinates is 
filled out according to the IDFILL* subroutine with the last element 
negative denoting the end of the string. The order specified by the 
vector must correspond to the coupling spring DOF's used to couple 
the two substructures together. 

KCPL “ Coupling spring stiffness matrix which connects only two substructures. 

The DOF's for the spring are in the order defined by ENTR2A and ENTR2B, 
respectively. * 9^). 


♦Explanation at end of section. 


Entry Point 3 
Card Number 
1 

$ENTR3 

NK 

IV 


Input Data 


NK^IV 


Format 

SENTR3 

Definitions 


NAMELIST Input for Card 1 

Nunrt)er of elements specified in "IV" vector 

Input vector defining how many modes are to be kept and reduced for 
all of the substructures. The elements are read in groups of three 
where the first element corresponds to soiro "ISUB", the second element 
to the number of modes to be "kept" for ISUB and the third element 
to the number of modes to be "reduced" for ISUB. Only the ISUB's 
specified by IV will be coupled together even though there may be 
several more substructures and springs residing in the basic data 
set. A maximum of 100 modes may be specified as "keot" and a maximum 
of 200 as "reduced". 


Entry Point 4 
Card Nun^er Input Data 

1 LAMDAI> 


Format 

$ENTR4 
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Definitions 


SENTR4 

LAMDA0 

Entry Point 5 
Card Number 
1 

$ENTR5 

NM0DES 

MDl 

MD2 

PFl 

PF2 

GMIN 

NTAPE 

KM((IDES 

LIST 


NAMELIST Input for Card 1 

p2 value used in the dynamic transformation for reduction. Where p is 
the circular frequency (Rad/Sec) and%(*p^) the eigenvalue. 


Input Data Format 
NMIjIDES. MDl. MD2, PF2, GMIN $ENTR5 
NTAPE. KM^OES. LIST 

Definitions 


NAMELIST Input for Card 1 

N{A.V.), the total number of physical eigenvectors to be calculated 
after coupling solution. N equals the total number of “kept" DOF. 

-1, do not calculate any physical modes. 

0 (A.V.), First mode number of eigenvectors to be printed. If MDl 
.EQ, 0, no modes are printed. 

N (A.V.), last mode number of eigenvectors to be printed. (MD2i.N) 

0 (A.V.), first participation factor mode number to be printed. If 
PFl, EQ. 0, no participation factors are printed. 

N (A.V.) last participation factor mode number to be printed. (PF2 ^N) 

0 (A.V.), minimum value for participation factors. Absolute values 
less than this will be set equal to zero and not printed. (Use WRITE 
Subroutine) . 

0 (A.V.), output DYNAMO tape containing coupled eigenvalues and eigen- 
vectors. If NTAPE ,EQ.O, no output tape is prepared. Tape will contain 
eigenvalues, frequencies in rad/sec, frequencies in Hz, and eigenvectors. 
Initialize tape, if necessary, on card 1, EP*0. 

N (A.V.), total number of modes and frequencies to be output on 
tape. (KMODESfN). 

0 (A.V.), LTAPE option ignored. 

1, a call to LTAPE (NTAPE) is initiated. 
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SUBROUTINE IDFILL (ID. IV, JV, NN. NKEPT) 


SPECIAL ROUTINE FOR GENERATING FULL DOF IDENTIFICATION VECTORS 
FOR INPUT TO SUBROUTINES REVSYM AND REV ADD. 

CALLS SUBROUTINE ZZBOMB 

PROGRAMMED BY EJ KUHAR NOVEMBER 1972 


SUBROUTINE ARGUMENTS 


c 

ID * 

INPUT VECTOR 

SIZE (M). WHERE M.LE.NN 

c 

IV = 

OUTPUT VECTOR EXPANDED FROM ID 

SIZE (NN) 

c 

JV = 

OUTPUT VECTOR ORDERED BY IV 

SIZE (NN) 

c 

NN = 

INPUT VECTOR LENGTH FOR IV , JV 




.GT.O FILLS FIRST NKEPT 

LOCATIONS OF IV FRW ID 


.LT.O FILLS LAST NN - NKEPT LOCATIONS OF IV FROM ID 
NKEPT = OUTPUT NUMBER OF ELEMENTS SPECIFIED BY ID, IF NN.LT.O 

NKEPT = ABS(NN) - NO. OF ELEMENTS SPECIFIED BY ID. 

C 

C**** ID VECTOR **** 

C THE ID VECTOR IS USED TO SPECIFY A PARTICULAR SEQUENCE OF INTEGERS 
C TO BE EXPANDED SEOUENTIALLY IN THE IV VECTOR. INCLUSIVE GROUPS OF 

C INTEGERS TO BE GENERATED IN ASCENDING ORDER MAY BE SPECIFIED AT ANY 

C ONE TIME BY USING THREE ELEMENTS OF ID WHERE 
C lOfI) = N1 

C ID (I+l) = 0 

C 10 (1+2) = N2 (N2, GT,N1+1) 

C THE INTEGERS FROM N1 TO N2 WILL BE SEOUENTIALLY LOADED INTO THE IV 

C MATRIX STARTING FROM THE FIRST AVAILABLE LOCATION IN IV. 

C IF THE TOTAL NUMBER OF INTEGERS SPECIFIED BY ID IS . LT . NN, THE LAST 

C ELEMENT OF ID MUST BE 'NEGATIVE'. NKEPT WILL BE CALCULATED AND THE 
C UNSPECIFIED INTEGERS WILL BE INSERTED INTO THE REMAINING IV LOCATIONS 
C 

C**** IV. J’ VECTORS **** 

C THE IV VECTOR WILL CONTAIN NN NON-ZERO ELEMENTS AND IS USUALLY 
C USED AS THE INPUT VECTOR 'IVEC IN SUBROUTINE REVSYM WHICH REVISES 
C THE DOF'S TO APPEAR IN THE ORDER SPECIFIED BY IVEC. 

C THE JV VECTOR REPRESENTS THE EQUIVALENT INPUT VECTOR FOR REVADD. 



ooooooooooooooooooooooooo ooooooooooo 


SUBROUTINE READ 


(A.NR,NC,KR,KC) 


READ MATRIX OF REAL NUMBERS FROM CARDS OR TAPE AND RPINT IT. WRITE 
MATRIX ON TAPE IF SO INDICATED (BY HAVING THE WRITE-TAPE NUMBER IN 
COLUMNS 79-80). 

THE EXPLANATION OF FORMATS USED BELOW IS . . . 

A - DENOTES ANY KEY PUNCH SYMBOL. (EG. M/K). 

I - DENOTES AN INTEGER NUMBER. (EG. 436). 

E - DENOTES A REAL NUMBER. (EG. 24.963). 

**** CARD INPUT **** 


FIRST CARD 


MATRIX NAME. NUMBER OF ROWS. NUMBER OF COLUMNS 
WITH A6.14.15 FOWAT. 

REMARKS IN COLUMNS 16-69. A-TYPE FORMAT. 

BLANK IN COL. 71 (NEW PAGE FOR TAPE WRITE PRINTOUT) 

* IN COL. 71 (SAME PAGE FOR TAPE WRITE PRINTOUT) 

$ IN COLUMN 72 FOR WRITE-TAPE INITIALIZATION. 
WRITE-TAPE CONTROL IN COLUMNS 73-78, MAY BE BLANK. OR 
THE VK)RDS REWIND OR LIST. OR (WHEN $ IN COLUMN 72) 

THE WRITE-TAPE-ID (EG, T1234). 

WRITE-TAPE NUMBER IN COLUMNS 79-80, (EG. 21). 


MIDDLE CARDS 


LAST CARD - 


DATA WITH FORMAT (215, 4E17) 

1- ST 15 IS THE ROW NUMBER, 

2- MD 15 IS THE COL. NUMBER OF THE NEXT El 7 FIELD. 
NEXT 4E17 ARE ELEMENTS OF THE MATRIX. 

TEN ZEROS IN COLUMNS 1-10. 


**** TAPE INPUT **** 


ONE CARD 


CALLS FORMA 


- MATRIX NAME, ZERO OR MINUS THE LOCATION NUMBER OF MATRIX 
ON READ-TAPE. READ-TAPE NUMBER (IF MINUS, NO PRINTOUT). 
MATRIX RUN NUMBER WITH A6,14,15,A6 FORMAT. 

- READ-TAPE CONTROL IN COLUMNS 22-27, MAY BE BLANK, OR THE 
WORDS REWIND OR LIST. 

- REMARKS IN COLUMNS 28-69, A-TYPE FORMAT. 

- $ IN COLUMN 72 FOR WRITE-TAPE INITIALIZATION. 

- WRITE-TAPE CONTROL IN COLUMNS 73-78, MAY BE BLANK, OR 
THE WORDS REWIND OR LIST, OR (WHEN $ IN COLUMN 72) 

THE WRITE-TAPE- ID (EG, T1234). 

- WRITE-TAPE NUMBER IN COLUMNS 79-80, (EG, 21). 

SUBROUTINES INTAPE, LTAPE. PAGEWD, RTAPE, WRITE .WTAPE.ZZBCMB 


C 

C 

C 

C 

C 

c 

c 


SUBROUTINE ARGUMENTS 

A = OUTPUT MATRIX READ FROM CARDS OR TAPE. 

NR = OUTPUT NUMBER OF ROWS IN MATRIX A. 

NC = OUTPUT NUMBER OF COLS IN MATRIX A. 

KR = INPUT ROW DIMENSION OF A IN CALLING PROGRAM. 

KC = INPUT COL DIMENSION OF A IN CALLING PROGRAM. 



4.2 OUTPUT DATA 


Output data is saved on a user-designated file/tape. It consists of coupled eigen- 
values, frequencies in rad/sec, frequencies in Hz, and the coupled eigenvector 
partitions with assigned matrix names of LAKSYS, IOMEGA, FREQ, PHYSl, PHYS2, . . . 
PHYSM, respectively. PHYSM is the last substructure partitioned eigenvector where 
Is the total number of substructure partitions. This data can be accessed via 
the “READ" subroutine. 

A data listing of the user tape can be obtained by using the "LTAPE" option in EP-5. 
4.3 DATA FILES 

Seven data files are required to execute SCAMP. These data files have been assigned 
the specific file codes of 11 through 17. The first five files (11-15) are used to 
save basic data. If the user desires to re-start SCAMP, these files must be saved 
on a data tape(s) after a SCAMP execution. For re-start, the user must load files 
11 through 15 with the old data. Files 16 and 17 are scratch files used by SCAMP. 
Additional files (or tapes) may be specified by the user for data input and/or output. 
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APPENDIX A 
SAMPLE PROBLEM 


This appendix contains a sample problem using all 5 entry 
points of SCAMP. The sample problem consists of a 20 degree of 
freedom longitudinal rod model shown in the sketch below. Two 

substructures are used, each with 10 degrees of freedom. The 

« 

output data printed by SCAMP follows a listing of the actual input 
data cards used to execute the program. 
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